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This  paper  develops  a  computational  model  to  represent  details  of  reactive  porous-media  transport,  ele¬ 
mentary  catalytic  chemistry,  and  electrochemistry  within  unit  cells  of  segmented-in-series  solid  oxide 
fuel  cell  (SIS-SOFC)  modules.  Because  the  composite  electrode  structures  are  thin  (order  of  tens  of 
microns),  electrochemical  charge-transfer  chemistry  can  proceed  throughout  the  composite  electrode 
structures.  Modeling  such  spatially  distributed  charge  transfer  is  significantly  more  complex  than  mod¬ 
eling  situations  where  the  charge  transfer  can  be  represented  at  an  interface  between  electrode  and 
electrolyte.  The  present  model  predicts  electric-potential  fields  of  electrode  and  electrolyte  phases,  with 
the  charge-transfer  rates  depending  upon  local  electric-potential  differences  and  the  local  gas-phase 
composition.  The  paper  summarizes  the  underpinning  physical  and  chemical  models  and  uses  examples 
to  illustrate  and  interpret  important  aspects  of  SIS  performance. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fig.  1  illustrates  the  layout  of  a  segmented-in-series  (SIS)  solid- 
oxide  fuel  cell  (SOFC)  module  in  which  planar  SIS  cells  are  arrayed 
as  narrow  strips  on  the  outside  of  a  porous  support  structure. 
Fuel  flows  inside  the  support  and  the  cathode  sides  of  the  SIS 
are  exposed  to  air  on  the  outside.  Each  cell  is  composed  of  a 
membrane-electrode  assembly  (MEA)  and  an  interconnect  struc¬ 
ture.  The  MEA  consists  of  a  cermet  anode  (e.g.,  Ni-yttria-stabilized 
zirconia,  Ni-YSZ),  dense  electrolyte  (e.g.,  YSZ),  and  composite  cath¬ 
ode  (e.g.,  strontium-doped  lanthanum  manganate-YSZ,  LSM-YSZ). 
The  cathode  in  the  present  study  is  a  two-layer  design  using  a  pure 
LSM  current  collection  layer  above  an  LSM-YSZ  functional  layer. 
Each  layer  is  usually  a  few  tens  of  microns  thick  and  the  support 
structure  is  on  the  order  of  a  millimeter  thick.  The  width  of  each  cell 
is  typically  on  the  order  of  a  few  millimeters.  As  illustrated  here,  a 
porous  catalyst  layer  at  the  interface  between  the  support  structure 
and  the  fuel  channel  can  promote  reforming  hydrocarbon  fuels. 

Oxygen  ions  O2-  are  produced  within  the  cathode  by  the  electro¬ 
chemical  reduction  of  gas-phase  oxygen,  using  electrons  supplied 
from  the  anode  of  the  neighboring  cell.  Electric  current  is  produced 
along  the  width  of  the  MEA  as  the  oxygen  ions  that  are  transported 
across  the  dense  electrolyte  are  electro-oxidized  by  the  fuel  H2. 
Oxygen  ions  are  transported  across  the  dense  electrolyte  within 
the  MEA,  but  the  cells  are  isolated  such  that  oxygen  ions  cannot 
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be  transported  between  cells.  The  inset  in  Fig.  1  shows  that  each 
cell  is  connected  electrically  in  series  with  the  lateral  neighboring 
cells.  The  electrons  produced  within  the  anode  of  one  cell  are  con¬ 
ducted  into  the  cathode  of  the  neighboring  cell.  Electric  current 
from  the  anode  interconnect  at  one  end  of  the  module  is  delivered 
to  an  external  circuit,  with  electrons  being  returned  to  the  cathode 
interconnect  at  the  other  end  of  the  module.  At  the  full  system  level, 
multiple  modules  can  be  connected  electrically  either  in  series  or 
parallel. 

1.1.  Potential  advantages  of  SIS 

Segmented-in-series  architectures  offer  several  potential  ben¬ 
efits.  Electrical  current  is  carried  laterally  through  width  of  each 
cell,  which  can  produce  significant  ohmic  losses,  especially  in 
the  cathode  where  the  materials  can  have  relatively  high  elec¬ 
trical  resistance  R.  However,  connecting  many  cells  electrically  in 
series  builds  up  module  voltage,  thus  reducing  ohmic  (/2R)  losses 
because  the  high  voltage  reduces  current  I  for  a  certain  electrical 
power.  Overall  cell  performance  benefits  by  reducing  lateral  elec¬ 
trical  resistance,  which  is  usually  accomplished  by  using  short  cell 
widths.  One  objective  of  the  present  paper  is  to  determine  optimal 
cell  widths. 

The  SIS  approach  confines  Ni  to  the  thin  anode  structure,  thus 
reducing  cost  by  significantly  reducing  the  amount  of  Ni  required  in 
comparison  to  more  common  Ni-YSZ  anode-supported  SOFC  archi¬ 
tectures.  The  SIS  supports  do  not  carry  electrical  or  ionic  current, 
and  thus  can  be  fabricated  from  relatively  low-cost  materials  such 
as  partially  stabilized  zirconia  (PSZ).  The  MEA  layers  can  be  fab- 
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Fig.  1.  Section  of  a  segmented-in-series  module. 


ricated  using  relatively  low-cost  screen-printing  technology.  The 
model  presented  in  this  paper  provides  quantitative  tool  to  assist 
evaluating  design  alternatives. 

1.2.  Prior  literature 

Alternative  SIS-SOFC  architectures  have  been  studied  since  the 
mid  1960s  [1-3].  Initial  efforts  considered  cells  that  were  applied 
as  rings  to  the  exterior  of  cylindrical  porous  support  tubes.  Later 
designs  considered  planar  cells  that  were  applied  to  flat  porous  sup¬ 
port  structures  (e.g.,  Fig.  1).  Mitsubishi  Fleavy  Industries  continues 
to  develop  SIS  architectures  on  circular  tubes  [4]  and  Rolls-Royce 
is  developing  planar  architectures  [5,6].  Rolls-Royce  uses  the  name 
Integrated  Planar  Solid  Oxide  Fuel  Cell  (IP-SOFC)  to  describe  their 
design.  Recently,  Liu  et  al.  have  fabricated  cone-shaped  anode- 
supported  SIS-SOFCs  with  good  performance  [7,8]. 

Barnett  et  al.  have  published  a  number  of  studies  concerning  sys¬ 
tem  design  and  optimization  [9-12].  Zhan  and  Barnett  suggested 
using  an  inner  reforming  or  partial-oxidation  catalyst  (as  illustrated 
in  Fig.  1 )  to  enable  the  use  of  hydrocarbon  fuels.  The  porous  sup¬ 
port  structure  also  serves  as  a  barrier  layer  that  facilitates  the  use  of 
hydrocarbon  fuels  without  carbon  deposits  [13,14].  Based  on  mod¬ 
eling  and  experimental  efforts,  Barnett  et  al.  have  achieved  high 
power  density  by  reducing  the  lateral  resistance  losses  across  the 
electrodes  and  maintaining  large  active  cell  area  through  lower¬ 
ing  the  cell  width  and  interconnect  area  [9,11,12].  Lai  and  Barnett 
[9]  presented  a  model  to  analyze  SIS  performance  as  functions  of 
cell  and  interconnect  geometry,  support  material,  cell  area-specific 
resistance,  and  electrode  sheet  resistance.  These  results  showed 
that  power  densities  around  1  W  cm-2  can  be  obtained  using  cell 
widths  of  about  1  -2  mm.  Lai  and  Barnett  also  reported  that  the  cath¬ 
ode  electrical  resistance  can  be  significantly  reduced  by  applying  a 
porous  LSM  current-collection  layer  above  the  cathode  [11].  Using 
Ni-YSZ|  YSZ|Pt-YSZ  cells  operating  on  humidified  H2  at  800 0  C  and 
atmospheric  pressure,  Pillai  et  al.  [12]  measured  power  densities 
of  approximately  0.7  W  cm-2  using  cell  widths  of  approximately 
1.2  mm  and  interconnect  widths  of  approximately  0.2  mm. 

Flaberman  and  Young  have  developed  three-dimensional 
computational-fluid-dynamics  (CFD)  models  to  investigate  the 
effects  of  fuel  and  air  flow  as  well  as  heat  and  mass  transport  at  the 
system  level  for  IP-SOFC  designs  [15-17].  Several  other  groups  have 
also  developed  planar  SIS  models  that  incorporate  porous-media 
transport,  coupled  with  reforming  chemistry  and  electrochemistry 
[18,19].  Cui  and  Cheng  developed  two-dimensional  axisymmetric 
models  for  the  tubular  SIS-SOFC  to  analyze  the  effects  of  the  cell 
geometry  on  the  electrical  performance  and  species  transport  [20]. 


Compared  to  prior  literature,  the  present  model  makes  signifi¬ 
cant  advances  in  the  fundamental  representation  of  chemistry  and 
electrochemistry.  Electric  potentials  for  both  ion-  and  electron¬ 
conducting  phases  are  modeled  throughout  the  entire  cell.  Lienee, 
both  ionic  and  electron  fluxes  are  predicted  throughout  the  system. 
Electrochemical  charge-transfer  chemistry  depends  on  the  local 
temperature,  gas-phase  composition,  and  electric-potential  differ¬ 
ences  between  phases.  The  spatial  extent  of  the  charge-transfer 
region  depends  on  electrode  structure,  including  primary  particle 
sizes,  phase  densities,  porosity,  tortuosity,  etc.  Reactive  porous- 
media  gas-phase  transport  is  represented  using  a  Dusty-Gas  model. 
The  model  also  represents  elementary  catalytic  chemistry  (typi¬ 
cally  tens  of  elementary  reactions)  within  the  composite  anode  and 
within  the  catalytic  layer  at  the  interface  between  the  fuel  flow 
and  the  support  structure.  The  elementary  catalytic  chemistry  rep¬ 
resents  the  internal  reforming  or  partial  oxidation  of  hydrocarbon 
fuels  as  well  as  other  thermal  chemistry  such  as  water-gas-shift 
process. 

2.  Physical  and  chemical  model 

This  paper  focuses  on  the  design  and  performance  of  individ¬ 
ual  unit  cells,  including  the  interconnect  region  and  the  underlying 
support  structure.  Flowever,  the  fuel  flow  within  the  channels  and 
the  exterior  air  flow  are  not  directly  modeled.  The  model  is  formu¬ 
lated  in  terms  of  continuum  differential  equations  that  describe  the 
electric  potentials  for  the  electrode  and  electrolyte  phases  as  well  as 
Faradaic  charge  transfer  between  phases.  The  porous-media  trans¬ 
port  of  gases  within  the  pore  spaces  is  modeled  with  a  Dusty-Gas 
model.  Catalytic  reforming  and  partial  oxidation  of  fuels  is  based 
upon  elementary  reaction  mechanisms. 

Fig.  2  illustrates  the  structure  of  the  unit  cell,  including  nomi¬ 
nal  dimensions.  The  figure  calls  attention  to  a  “unit  width”  and  an 
“active  width”.  The  unit  width  is  the  width  of  the  cell  repeating 
unit,  which  includes  interconnect  and  cell-isolation  features.  The 
active  width  is  the  width  over  which  the  anode  and  cathode  struc¬ 
tures  overlap.  Electrical  current  is  produced  primarily  within  the 
active  width.  As  illustrated  in  Fig.  2,  the  present  paper  considers 
the  lateral  widths  of  features  within  the  interconnect  regions  to  be 
fixed  at  80  p,m.  Assuming  a  fixed  interconnect  structure,  the  model 
is  used  to  determine  unit  width  that  maximizes  cell  performance. 

Because  the  anode  and  cathode  thicknesses  are  on  the  order  of 
50  |xm,  the  electrochemical  charge-transfer  processes  are  likely  to 
be  distributed  throughout  most  of  the  porous  electrode  structure 
[21].  Therefore,  SIS  models  should  accommodate  charge-transfer 
electrochemistry  throughout  the  MEA.  Because  a  detailed  discus- 
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Fig.  2.  Illustration  of  the  MEA  structure  of  the  SIS  unit  cell,  showing  nominal  physical  dimensions. 


sion  of  the  distributed  charge  transfer  model  is  provided  elsewhere 
[21  ],  the  transport  and  chemistry  theory  is  presented  here  in  sum¬ 
mary  form. 


2.1.  Porous-media  transport  and  chemistry 


Chemically  reacting  gas-phase  species  transport  through  the 
porous  electrodes  is  represented  by  species  and  overall  mass- 
continuity  equations  as 


d(4>gPYk) 

dt 

d(4>sp) 

dt 


+  v  •  jk  =  wksk, 

(i) 

Kg  Kg 

k=  1  k=  1 

(2) 

The  independent  variables  are  time  t  and  the  spatial  coordinates. 
The  dependent  variables  are  the  gas-phase  mass  density  p  and 
the  mass  fractions  Yk.  Other  parameters  include  the  porosity  </>g 
and  species  molecular  weights  Wk.  Gas-phase  species  mass  fluxes 
j/<  through  the  porous  structure  are  determined  from  the  Dusty- 
Gas  model  (DGM)  [22,23],  which  is  an  implicit  relationship  among 
the  gas-phase  species  molar  fluxes  J/0  molar  concentrations  [Xk], 
concentration  gradients,  and  the  pressure  gradient  as 


WXiUfc-M  ,  h 

ffk  [*T  He  DeKn 


-V[Xk] 


[xk]  fh 

Dt,Kn  * 


vp, 


(3) 


where  /x  is  the  mixture  viscosity,  [XT]  =p/RT  is  the  total  molar  con¬ 
centration,  p  is  pressure,  and  Bg  the  permeability.  The  mass  fluxes 
j/<  are  related  to  the  molar  fluxes  ]k  as  =  WjJk.  This  formulation 
considers  molecular  transport  via  ordinary  and  Knudsen  diffusion 
as  well  as  pressure-driven  Darcy  flow.  The  effective  binary  and 
Knudsen  diffusion  coefficients  can  be  evaluated  as 


Du  = 


kit 


_  2  rp0g 
k*n  ~  3  tg 


8  RT 


(4) 


which  are  related  to  the  ordinary  multicomponent  diffusion  coeffi¬ 
cients  Dki  via  the  porous-media  microstructure,  including  average 
pore  radius  rp  and  tortuosity  rg.  For  the  problems  considered  in  the 
present  paper,  pressure  gradients  within  the  porous  structures  are 
very  small,  making  the  pressure-driven  species  transport  negligi¬ 
ble.  Nevertheless,  the  full  Dusty-Gas  model  is  retained. 

The  molar  production  rates  of  gas-phase  species  due  to 
heterogeneous  thermal  catalytic  chemistry  and  electrochemical 
charge-transfer  reactions  are  represented  by  sk.  Because  the  tem¬ 
peratures  are  relatively  low  (^800  °C)  and  pore  size  within  an 
electrode  is  comparable  to  the  molecular  mean-free-path  length, 
homogeneous  gas-phase  kinetics  can  be  neglected.  The  species 
production  rates  sk  are  functions  of  temperature,  gas  composi¬ 
tion,  surface-species  coverages,  and  electric-potential  differences 


between  the  electrode  and  electrolyte  phases.  In  addition  to  gas- 
phase  species,  the  chemistry  also  depends  upon  surface-adsorbed 
species.  The  temporal  variations  of  site  coverages  0k  of  surface- 
adsorbed  species  are  represented  as 

=  %  (k  =  1, . . . ,  Ks),  (5) 

where  r  is  available  site  density  and  I<s  is  the  number  of  the  surface- 
adsorbed  species.  At  steady  state  the  net  production  rates  of  surface 
species  vanish  (i.e.,  sk  =  0). 

Heterogeneous  reforming  within  the  Ni-YSZ  anode  structure 
represented  by  an  elementary  reaction  mechanism  that  incor¬ 
porates  steam  and  dry  reforming  as  well  as  partial  oxidation 
on  the  Ni  surfaces  [24,25].  The  reaction  mechanism  also  repre¬ 
sents  the  important  effects  of  water-gas-shift  processes  (globally, 
H20  +  C0^H2  +  C02).  In  addition  to  reforming  within  the  Ni-YSZ 
anode,  a  Rh  catalyst  layer  at  the  interface  between  the  support 
structure  and  the  fuel  channel  may  be  used  to  promote  hydrocarbon 
reforming.  The  catalytic  reforming  of  CH4  on  Rh  is  also  represented 
by  an  elementary  heterogeneous  reaction  mechanism  [26]. 

2.2.  Charge  conservation 

Charge  transport  and  charge  transfer  both  depend  upon  spa¬ 
tial  electric-potential  variations  within  the  electron-conducting 
phase  (i.e.,  the  Ni  phase  in  the  anode  and  LSM  in  the  cathode) 
and  ion-conducting  phase  (i.e.,  the  ceramic  YSZ  in  the  anode,  cath¬ 
ode  and  electrolyte).  Because  the  composite  electrodes  are  thin, 
quantitative  predictions  require  modeling  electric-potential  vari¬ 
ations  throughout  the  entire  MEA.  This  is  in  contrast  to  modeling 
anode-supported  cells,  where  it  is  sometimes  appropriate  to  model 
the  active  charge-transfer  region  without  spatial  resolution  at  the 
interface  between  the  composite  anode  and  the  dense  electrolyte. 
In  a  Ni-YSZ|  YSZ|LSM-YSZ  MEA,  three  electric-potential  fields  are 
relevant.  They  are  the  electric  potential  for  the  electron-conducting 
phase  within  the  composite  anode  @a,  the  electric  potential  for  the 
electron-conducting  phase  within  the  cathode  @c,  and  the  elec¬ 
tric  potential  for  the  electrolyte  phase  which  spans  the  anode, 
dense  electrolyte,  and  cathode.  The  three  electric  potentials  are 
governed  by  charge-conservation  equations.  For  steady-state  oper¬ 
ating  conditions  the  charge-conservation  equations  may  be  stated 
as [21] 


'  sa,e  within  anode 


V.(creeV<Z>e)=  < 

0 

within  electrolyte , 

(6) 

,  Sc,e 

within  cathode 

V(<7fV<Z>a)  =  - 

Sa,e 

within  anode, 

(7) 

V  ■  (erfV^c)  =  - 

Sc,e 

within  cathode. 

(8) 
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In  these  equations  erf  and  cr£,  are  the  effective  conductivities  of 
the  electron-conducting  phases  within  the  composite  anode  (i.e. 
Ni)  and  cathode  (i.e.,  LSM),  respectively.  The  effective  ion  conduc¬ 
tivity  for  the  ion-conducting  phase  (i.e.,  YSZ)  is  a £.  The  effective 
conductivities  of  the  composite  electrodes  can  be  significantly  dif¬ 
ferent  from  the  pure  material  properties.  Generally  speaking  the 
ion  conductor,  the  electron  conductor,  and  the  pores  each  occupy 
approximately  one  third  of  the  composite  structure.  Thus,  on  a  sim¬ 
ple  geometric  basis,  one  expects  that  the  effective  conductivity  to 
be  no  greater  than  one  third  of  the  pure  material  conductivity. 
Because  of  necking  between  particles,  particle-particle  interface 
resistances,  tortuous  conduction  paths,  and  material  inhomo¬ 
geneities  such  as  grain  boundaries,  the  effective  conductivities  are 
usually  much  lower  than  might  be  expected  by  considering  the 
phase  volume  alone.  Effective  properties  can  be  estimated  using 
percolation  theory  [27]  and  derived  from  detailed  simulations  of 
synthesized  microstructures  [28]. 

Electrochemical  transfer  of  charge  between  phases  affects  the 
electric  potentials  within  the  phases.  In  the  continuum  conser¬ 
vation  equations  the  net  charge-transfer  rates  are  represented  as 
source  terms  on  the  right-hand  sides  of  Eqs.  (6)-(8).  The  nomencla¬ 
ture  sm,e  has  the  meaning  that  the  phase  m  may  be  either  the  anode 
“a”  or  the  cathode  “c.”  Charge-transfer  reactions  usually  involve 
both  charged  and  uncharged  species.  The  second  subscript  “e”  in 
sm>e  indicates  the  transfer  of  electrons,  which  are  the  only  rates 
that  are  directly  relevant  to  the  charge-conservation  equations. 
A  positive  value  of  sm,e  indicates  the  transfer  of  negative  charge 
into  phase  m,  which  tends  to  decrease  the  electric  potential  of  the 
phase.  Note  that  the  charge-transfer  source  terms  appear  in  the 
charge-conservation  equations  for  both  participating  phases,  but 
with  opposite  signs.  That  is,  all  the  charge  that  leaves  one  phase 
enters  the  partner  phase.  Further  details  about  the  charge-transfer 
formulation  is  found  in  Zhu  and  Kee  [21  ]. 

2.3.  Butler-Volmer  kinetics 


equilibrium  electric  potential,  the  chemistry  proceeds  at  equal 
and  opposite  rates  in  the  anodic  and  cathodic  directions.  The 
Butler-Volmer  expression  describes  the  net  rate  of  a  reversible 
reaction.  The  exchange  current  density  represents  the  magnitude 
of  the  equal  and  opposite  current  densities  at  equilibrium  elec¬ 
tric  potential.  The  value  of  z0  depends  upon  the  length  of  three 
phase  boundary,  the  effectiveness  of  electrocatalyst,  the  temper¬ 
ature,  and  the  activities  of  the  participating  species.  Because  the 
Butler-Volmer  formulation  represents  a  global  charge-transfer 
process,  the  symmetry  factors  a  are  not  the  same  as  they  would 
be  for  an  elementary  reaction.  Moreover,  aa+ac  ^  1  as  would  be 
the  case  for  an  elementary  reaction. 

The  activation  overpotential  rjact  is  defined  as 

Pact  =£a- £aeq  =  (*a  -  <*>e)  -  (<Z>a  -  <Z>e)eq.  (12) 


That  is,  Ea  =  {@a-  is  the  electric-potential  difference  between 
the  electrode  phase  and  the  electrolyte  phase,  and  £fq  = 
{<Pa-0e)eq  is  the  electric-potential  difference  that  causes  the 
charge-transfer  reaction  to  be  equilibrated  (i.e.,  proceeding  in  the 
anodic  and  cathodic  directions  at  equal  and  opposite  rates). 

Zhu  et  al.  [23,21,32]  developed  expressions  for  the  anode  and 
cathode  exchange  current  densities  as 


1o,h2  =  i'h2 


(PH2/pt,J(1““a/2)(PH20)“a/2 


KH2 

1  +(Ph2/P&2) 


1/2 


(13) 


'o,o2  =  <02 


(po2/Po2ra/2 

1  +  (Po2/Pq2)1/2  ' 


(14) 


where  the  gas-phase  partial  pressures  pk  are  measured  in  atmo¬ 
spheres.  Zhu,  et  al.  derive  expressions  for  the  factors  p^2  and 
Pl j2,  which  depend  upon  heterogeneous  adsorption  and  desorption 
rates  [23]. 

Temperature  dependencies  for  the  exchange  current  densities 
are  expressed  as 


Within  a  composite  electrode  (e.g.,  Ni-YSZ  anode  or  LSM-YSZ 
cathode)  the  overall  charge-transfer  processes  involve  a  com¬ 
plex  interaction  of  thermal  chemistry  and  electrochemistry.  Even 
hydrogen  electro-oxidation  within  a  Ni-YSZ  anode  can  involve  a 
complex  set  of  elementary  reactions  [29-33].  The  present  paper 
represents  charge-transfer  kinetics  in  a  Butler-Volmer  form  that  is 
derived  from  elementary  reactions  [23,32]. 

Within  the  Ni-YSZ  anode,  H2  is  assumed  to  be  the  only  elec- 
trochemically  active  fuel.  Although  the  fuel  stream  may  contain 
CO  and  methane,  the  catalytic  reforming  and  water-gas-shift  pro¬ 
cesses  are  dominant  in  the  production  of  H2.  The  electrochemical 
oxidation  of  H2  within  the  Ni-YSZ  anode  and  02  reduction  within 
the  LSM-YSZ  cathode  can  be  represented  globally  as 

H2(g)  +  02-(el)^H20(g)  +  2e-(a),  (9) 

(l/2)02(g)  +  2e-(c)^02-(el),  (10) 

where  02-(el)  are  oxygen  ions  within  the  ionic-conducting 
electrolyte  phase,  e-(a)  are  the  electrons  within  the  electronic- 
conducting  anode  phase,  and  e-(c)  are  the  electrons  within  the 
electronic-conducting  cathode  phase.  In  Butler-Volmer  form  the 
charge-transfer  rates  (i.e.,  current  density  zBv,e)  can  be  expressed  in 
terms  of  the  local  activation  overpotential  pact  as 


*BV,e  =  *0 


exp 


/  O'aEpact  \ 

\~~rT~  ) 


-  exp 


/  oicFrjact  \ 

V  kT~  ) 


(ii) 


where  z0  is  the  exchange  current  density,  aa  and  ac  are  anodic 
and  cathodic  symmetry  factors.  The  first  term  in  the  square  brack¬ 
ets  represents  the  anodic  current  density  (producing  electrons) 
and  the  second  term  represents  the  cathodic  current  density.  At 
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where  Eak  is  activation  energies,  and  the  parameter  i*ef  k  is  assigned 
empirically  to  fit  measured  polarization  data  at  the  reference  tem¬ 
perature  Tre  f. 

The  spatially  varying  local  electric-potential  differences  within 
the  porous  anode  and  cathode  structures  are  defined  as 


Ea  =  0 a-<Z>e,a,  £c  =  &c  ~  #e.c,  (17) 

where  0e,m  represents  the  electrolyte-phase  electric  potential  in 
the  composite  anode  “a”  or  cathode  “c”.  The  local  electric  potential 
fields  are  obtained  directly  from  the  solutions  to  Eqs.  (6)-(8).  The 
activation  overpotentials  in  Eq.  (11)  for  the  anode  and  cathode  are 
defined  as 


Pact, a  —  Ea  —  Eaq,  Pact.c  —  Ec  —  Ecq, 


(18) 


which  depend  upon  the  local  equilibrium  electric-potential  differ¬ 
ences  £fq  and  £j:q.  Assuming  bulk  O2-  concentration  within  the  YSZ 
lattice  is  spatially  uniform  throughout  the  MEA  structure,  £fq  and 
£^q  can  be  evaluated  as  [21  ] 


Feq  _  ^H20  ^H2  RT  PH20,a 

3  “  2F  +  2 F  ln  PH2,a 


(19) 


^  RT 
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where  /x°  are  standard-state  chemical  potentials  and  pk  are  the 
local  gas-phase  partial  pressures  (in  atmospheres). 

2.4.  Boundary  conditions 

Boundary  and  interface  conditions  are  required  to  solve  the 
coupled  system  of  Eqs.  (1),  (2)  and  (6)-(8).  Fig.  3  illustrates  the 
numerous  boundary  and  interface  conditions  that  must  be  imposed 
throughout  the  multicomponent  cell  structure. 

At  the  interface  with  the  fuel  and  air  compartments  the  gas- 
phase  composition  is  assumed  to  be  that  within  either  the  bulk 
fuel  or  air  flow.  At  the  interface  between  the  cathode  and  the  air 
compartment  the  ionic  fluxes  vanish  (i.e.,  n  •  o^V<2>e  =  0)  and  the 
cathode  electric  potential  @c  can  be  set  to  be  a  reference  electric 
potential  as  @°.  At  the  interface  between  the  anode  and  the  fuel 
compartment  the  ionic  fluxes  vanish  (i.e.,  n  •  cr%¥(Pe  =  0),  and  the 
anode  electric  potential  can  be  specified  as  @°a.  The  operating 
cell  potential  is  £cen  =  <P°C  -  0°.  The  electron  flux  at  the  interface 
between  the  composite  electrode  structures  and  the  dense  elec¬ 
trolyte  must  vanish  (i.e.,  n  •  or|V0a  =  0,andn  •  aceV0c  =  0)  because 
the  dense  electrolyte  membrane  is  assumed  to  be  a  purely  ionic 
conductor. 

3.  Numerical  algorithm 

The  governing  conservation  equations  form  a  coupled  system  of 
nonlinear  partial-differential  equations.  The  independent  variables 
are  the  spatial  coordinates  and  the  dependent  variables  are  the 
gas-phase  density  pg,  mass  fractions  Yk  (/<  =  1, . . .,  I<g),  surface  cov¬ 
erages  0k  (/<  =  1, . . .,  Ks),  electric  potentials  @m  within  the  electrode 
and  electrolyte  phases.  The  entire  cell  and  support  structure  is  dis¬ 
cretized  using  a  finite-volume  representation  on  a  cartesian  mesh. 
Component-specific  material  properties  and  reaction  chemistries 
are  used  within  each  of  the  components. 

An  operator-decomposition  technique  is  used  to  decouple  the 
gas-phase  transport  and  reaction  processes.  The  numerical  solution 
is  obtained  by  solving  sequentially  and  iteratively  three  submodels: 
(1)  the  electric-potential  distributions,  (2)  the  gas-phase  species 
transport  within  the  porous  electrodes,  and  (3)  the  chemical  and 
electrochemical  reactions.  The  resulting  system  of  nonlinear  alge¬ 
braic  equations  is  solved  with  a  modified  Newton  iteration,  with 
a  preconditioned  Krylov  subspace  iteration  used  to  solve  the  lin¬ 
earized  algebraic  equations  [34]. 

The  chemical  kinetics,  which  does  not  involve  spatial  opera¬ 
tors,  is  evaluated  within  each  finite  volume  by  a  time-relaxation 
algorithm,  using  Limex  to  solve  the  differential-algebraic  equations 
[35  ].  Thermodynamic  properties,  transport  properties,  and  reaction 
chemistry  are  managed  using  the  Chemkin  software  [36]. 

4.  Results  and  discussion 

The  present  study  considers  an  SIS  architecture  with  specific 
component  thicknesses  and  interconnect  structures,  but  with  dif¬ 
ferent  unit-cell  widths.  As  introduced  in  Fig.  2  physical  dimensions 
are  as  follows.  The  composite  cathode  consists  of  two  layers  -  a 
30  [xm  LSM-YSZ  functional  layer  adjacent  to  the  dense  electrolyte 
layer  and  a  30|xm  LSM  layer  to  increase  the  lateral  electronic 
conduction.  The  dense  YSZ  electrolyte  membrane  is  20  |xm  thick 
and  the  composite  Ni-YSZ  anode  is  50  |xm  thick.  The  MEA  struc¬ 
ture  rests  on  a  support  layer,  that  is  composed  of  450  pum  of 
a  porous  chemically  inert  structure  above  a  Rh  catalyst  layer 
that  is  50|xm  thick.  The  electronically  conducting  interconnect 
is  an  L-shaped  structure  that  consists  of  two  80  [xm  wide  seg¬ 
ments.  Because  the  interconnect  must  survive  in  both  oxidizing  and 
reducing  environments,  materials  selection  is  limited.  One  viable 


Table  1 

Nominal  physical  parameters. 


Parameters 

Value 

Units 

Anode,  Ni-YSZ 

Thickness  (La) 

50 

pan 

Porosity  (0g) 

0.35 

Ni  volume  fraction  (0Ni) 

0.23 

YSZ  volume  fraction  (0 Ysz) 

0.42 

Tortuosity  (rg) 

4.5 

Ni  particle  radius  (rNi) 

0.5 

pm 

YSZ  particle  radius  (rY sz) 

0.5 

pm 

Specific  catalyst  area  (As) 

1000 

cm-1 

Exchange  current  factor  (i*ef  ) 

48,000 

Acnrr3 

Activation  energy  (Ea,H2 ) 

120 

kj  mol-1 

Reference  temperature  (Tref) 

800 

°C 

Anodic  symmetry  factor  (aa) 

1.5 

Cathodic  symmetry  factor  (ac) 

0.5 

Cathode,  LSM-YSZ 

Thickness  (Lc) 

60 

pm 

Porosity  (0g) 

0.35 

LSM  volume  fraction  (<2>lsm) 

0.31 

YSZ  volume  fraction  (<2>ysz) 

0.34 

Tortuosity  (rg) 

4.00 

LSM  particle  radius  (tlsm) 

0.625 

pm 

YSZ  particle  radius  (rY sz) 

0.625 

pm 

Exchange  current  factor  (i*ef  ) 

7200 

Acm~3 

Activation  energy  (Ea, o2 ) 

130 

kj  mol-1 

Reference  temperature  (Tre f) 

800 

°C 

Anodic  symmetry  factor  (aa) 

0.55 

Cathodic  symmetry  factor  (ac) 

0.45 

Electrolyte,  YSZ 

Thickness  (Le i) 

20 

pm 

Anode  support  diffusion  layer 

Thickness  (Lad) 

450 

pm 

Porosity  (0g) 

0.35 

Tortuosity  (rg) 

4 

Solid  volume  fraction  (0S) 

0.65 

Solid  particle  radius  (rs) 

0.75 

pm 

Interconnect  (LSC) 

Electrical  conductivity  (crLSC) 

1 

Scm-1 

Anode  support  Rh  catalyst  layer 

Thickness  (Lac) 

50 

pm 

Porosity  (0g) 

0.35 

Tortuosity  (rg) 

4 

Solid  volume  fraction  (0S) 

0.65 

Solid  particle  radius  (rs) 

0.75 

pm 

Specific  catalyst  area  (As) 

1000 

cm-1 

choice  is  strontium-doped  lanthanum  chromite  (LSC),  but  it  has 
low  electronic  conductivity  (tr  <  1  S  cm-1 ).  The  lateral  isolation  gap 
between  cathodes  is  160|xm  wide,  spanning  80|xm  of  intercon¬ 
nect  and  80  |xm  of  dense  electrolyte.  Reducing  the  widths  of  these 
80  |xm  spaces  would  be  functionally  beneficial.  However,  dimen¬ 
sions  smaller  than  80  |xm  may  be  difficult  to  fabricate  reliably  with 
technology  such  as  screen-printing. 

The  nominal  unit-cell  width  in  the  present  study  is  taken  to 
be  1900  |xm,  with  the  active  width  (i.e.,  width  over  which  the 
cathode  and  anode  are  both  in  direct  contact  with  the  dense  elec¬ 
trolyte)  being  1500  |xm.  A  particular  objective  of  the  paper  is  to 
evaluate  cell  performance  as  a  function  of  varying  unit-cell  active 
width. 

Table  1  summarizes  the  physical  dimensions  and  provides 
numerous  other  physical  parameters  and  properties.  In  all  cases 
the  cathode  is  exposed  to  air  at  atmospheric  pressure.  However, 
two  fuel  streams  are  considered:  (1)  humidified  H2  (3%  H20)  and 
(2)  a  mixture  resulting  from  partial  steam  reforming  of  CH4.  The 
cell  is  assumed  to  be  operated  isothermally  800  °C.  Two  cell  oper¬ 
ating  potentials  are  considered:  £cen  =  0.50  V  and  Fcell  =  0.75  V.  The 
cell  potential  Eceli  is  defined  to  be  the  difference  between  the  cath¬ 
ode  interconnect  potential  (left  side  of  unit  cell)  and  the  anode 
interconnect  (right  side  of  the  unit  cell). 
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Yk  =  Air 
n*V$c  =  0 
n*V$P  =  0  I 


[<&c  =  $c° 

n-V<Fe  =  0 

n-Vn  -  0 


n*V<&a  =  0 
n-VVfc  =  0 


n*V$e  =  0 

n-VYfe  =  0 


n-V$e  —  0 

Yk  =  Fuel ' 

ruei 

n-V$e  =  0 

n-VYjfc  =  0 

n-V<E>a  =  0 

n-VFfc  =  0 
\ _ _ _ / 

n-V$e  =  0 

^  J 

Cell  potential:  Ece n  =  Gas-phase  mass  fraction:  Yk 

Cathode  phase  electric  potential:  <hc  Anode  phase  electric  potential:  d>a 

Electrolyte  phase  electric  potential:  d>e 

Fig.  3.  Boundary  and  interface  conditions  throughout  the  MEA  unit  cell. 


4.1.  Characterization  of  a  unit  cell 

Consider  first  a  fuel  that  is  the  result  of  equilibrium  reforma¬ 
tion  of  an  initial  mixture  of  60%  CH4  and  40%  H20  at  800°  C  and 
atmospheric  pressure.  Specifically,  the  fuel  is  a  mixture  of  66.1% 
H2,  21.8%  CO,  11.6%  CH4,  0.3%  H20,  and  0.2%  C02. 

It  may  be  noted  that  the  gas  composition  flowing  into  the  SOFC 
fuel  channel  could  be  such  that  carbon  deposits  would  be  formed, 
especially  on  Ni.  However,  the  SIS  architecture  (Fig.  1)  provides 
certain  advantages  that  counteract  the  propensity  for  coking.  The 
chemically  inert  porous  support  structure  acts  as  a  barrier  layer 
[13,14]  that  impedes  the  transport  of  electrochemically  generated 
H20  from  the  Ni-YSZ  anode  and  impedes  the  transport  of  hydro¬ 
carbons  toward  the  anode.  Thus,  the  gas  mixture  within  the  Ni-YSZ 
anode  is  less  likely  to  coke  than  would  be  the  gas  composition  the 
flow  channel. 

Using  the  nominal  cell  dimensions  and  material  properties 
(Table  1 ),  Fig.  4  shows  predicted  gas-phase  mole-fraction  (percent) 
contours  within  the  pore  volumes  of  the  anode  and  support  layers. 
Qualitatively,  the  concentrations  of  fuel  species  (i.e.,  H2,  CO,  and 
CH4)  decrease  from  their  mole  fractions  within  the  fuel  channel 
toward  the  anode  structure  where  they  are  consumed.  The  concen¬ 
trations  of  product  species  (i.e.,  H20  and  C02 ),  which  are  generated 
within  the  anode,  decrease  as  they  are  transported  toward  the  fuel 
channel.  Concentration  gradients  are  much  higher  within  the  cat- 
alytically  active  anode  layer  than  they  are  in  the  chemically  inert 
support  structure.  Because  the  fuel  mixture  is  already  at  equilib¬ 
rium  composition  and  the  Rh  catalyst  layer  is  only  50  pm  thick,  the 
catalyst  layer  has  little  influence. 

Over  much  of  the  anode-layer  width  the  species  contours 
remain  nearly  parallel.  In  other  words,  the  chemistry  and  transport 
within  much  of  the  MEA  is  nearly  one-dimensional  in  the  direction 
normal  to  the  dense  electrolyte  layer.  However,  near  the  inter¬ 
connect  region,  where  electrochemistry  is  not  active,  significant 
two-dimensional  behavior  is  evident.  Geometric  design  and  mate¬ 
rials  selection  within  the  interconnect  region  is  one  of  the  more 
challenging  aspects  of  SIS  development. 

Electrochemical  variations  in  the  interconnect  region  affect  the 
species  fields  within  the  support  structure.  By  assumption,  the 


species  fluxes  within  the  support  structure  vanish  at  the  interfaces 
between  the  upstream  and  downstream  unit  cells.  Because  fuel  is 
consumed  along  the  length  of  the  fuel  channels,  this  assumption  is 
not  strictly  valid.  However,  the  zero-flux  assumption  at  the  unit¬ 
cell  boundaries  is  a  reasonable  approximation.  With  small  unit-cell 
widths  (e.g.,  1900  |jim),  relatively  little  electrical  power  is  gener¬ 
ated  within  each  unit  cell.  Consequently,  relatively  little  fuel  is 
consumed  and  products  generated  within  each  unit  cell.  Within 
a  unit  cell,  the  largest  species  fluxes  between  the  MEA  and  the  fuel 
channel  are  near  the  center  of  the  support  structure. 

In  addition  to  species  contours,  Fig.  4  also  shows  electric- 
potential  contours  within  the  electrolyte  phase  (YSZ)  and  electron 
pathlines  within  the  electron-conducting  phases  (Ni  and  LSM) 
throughout  the  MEA  structure.  As  with  the  species  contours,  the 
electric-potential  contours  are  nearly  one-dimensional  over  much 
of  the  MEA  width.  Most  of  the  electric-potential  variation  is  across 
the  thickness  of  the  dense-electrolyte  layer.  The  concentration  of 
electron  pathlines  in  the  LSM  current-collection  layer  is  the  result  of 
higher  electrical  conductivity  in  that  layer  compared  to  the  under¬ 
lying  LSM-YSZ  active  cathode.  Improving  lateral  conduction  within 
the  cathode  structure  is  an  important  consideration  in  improving 
overall  cell  performance.  Note  that  the  electron  pathlines  do  not 
cross  the  dense  electrolyte,  which  is  a  pure  oxygen-ion  conductor. 
Rather,  electrons  are  consumed  by  oxygen  electro-reduction  on  the 
cathode  side  and  produced  by  electro-oxidation  of  the  fuel  on  the 
anode  side. 

4.1.1.  Transverse  profiles 

Fig.  5  shows  profiles  of  electric  potentials  and  current  densi¬ 
ties  through  the  MEA  structure  at  three  lateral  positions  (165  pm, 
900  |jim,  and  1500  pan).  In  all  cases  the  cell  is  operating  at 
Eceii  =  0.75  V.  The  graphs  are  oriented  such  that  the  anode  is  on  the 
left  and  the  cathode  is  on  the  right.  Consider  first  Fig.  5a,  showing 
electric  potential  profiles.  The  solid  lines  represent  the  electrode 
phases  (Ni  and  LSM)  and  the  dashed  lines  represent  the  electrolyte 
phase  (YSZ).  The  electric-potential  profiles  are  nearly  overlapping, 
indicating  approximately  one-dimensional  behavior.  The  electrode 
electric  potential  is  nearly  uniform  within  the  Ni  at  <2>a  «  -  0.12  V 
and  within  the  LSM  at  0.78  V.  With  Ecen  =  0.75  V,  the  anode 
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Fig.  4.  Contour  maps  of  gas-phase  species  mole  fractions  within  the  anode  and  support  layers.  Electric-potential  contours  within  the  electrolyte  (YSZ)  phase  through  the 
MEA  structure.  Electron  path  lines  through  the  anode,  cathode  and  interconnect  structures.  To  assist  visualization  the  vertical  direction  is  expanded  by  a  factor  of  two. 


interconnect  is  fixed  at  <2>a  =  0  V  and  the  cathode  interconnector  is  at 
@c  =  0.75  V.  The  fact  that  the  electrode  electric-potential  difference 
is  greater  than  Ecell  (i.e.,  the  difference  in  interconnect  potentials) 
is  the  result  of  the  electric-potential  gradients  that  are  required  to 
support  lateral  current  flux. 

Within  the  LSM-YSZ  cathode,  the  electric  potentials  of  both 
the  cathode  phase  @c  and  electrolyte  phase  @e  increase  from  the 
cathode-air  interface  toward  the  interface  between  the  cathode 
and  the  dense  electrolyte.  Both  the  negatively  charged  ions  and 
electrons  are  being  transported  from  the  cathode-air  interface 
toward  the  interface  between  the  cathode  and  the  dense  elec¬ 
trolyte.  Such  charged-species  transport  proceeds  in  the  direction 
of  the  negative  electrochemical-potential  gradient.  In  other  words, 
the  negatively  charged  species  generally  are  transported  up  the 
electric-potential  gradient  (i.e.,  from  regions  of  relatively  more  neg¬ 
ative  electric  potential  toward  regions  of  relatively  more  positive 
electric  potential).  The  electric  potential  of  the  electrode  phase 
(LSM)  is  always  positive  relative  to  the  electric  potential  of  the 
electrolyte  phase  (YSZ).  That  is  local  double  layers  between  elec¬ 
trode  and  electrolyte  phases  are  such  that  A<2>  =  <2>LSm  -  ^ysz  >0. 


The  magnitude  of  the  double  layer  (i.e.,  A@)  is  less  positive  near 
the  dense-electrolyte  interface  than  it  is  deeper  into  the  cathode 
structure.  Consequently,  the  cathodic  charge-transfer  rate  is  high¬ 
est  near  the  dense-electrolyte  interface.  As  the  magnitude  of  A@ 
decreases  (i.e.,  the  YSZ  particles  become  relatively  more  positive) 
the  rate  of  the  charge-transfer  reaction  (globally  (l/2)02  (g)  +  2 
e-  (LSM)^O2-  (YSZ))  increases  (Fig.  5b).  This  can  be  understood 
qualitatively  in  the  sense  that  as  the  strength  of  the  double  layer 
decreases,  it  becomes  easier  to  transfer  negative  charge  from  the 
electrode  into  the  relatively  more  negative  electrolyte.  Within  the 
cathode  structure,  the  activation  overpotential  is  always  negative. 
In  the  context  of  the  Butler-Volmer  equation,  a  negative  overpo¬ 
tential  drives  the  charge-transfer  reaction  in  the  cathodic  direction 
(i.e.,  consuming  electrons). 

The  electric-potential  profiles  within  the  composite  anode 
increase  slightly  from  the  interface  between  the  anode  and  the 
dense  electrolyte  interface  toward  the  interface  between  the  anode 
and  support  structure.  As  in  the  cathode,  negatively  charged  ions 
and  electrons  are  transported  up  the  electric-potential  gradients. 
However,  because  the  Ni  electrical  conductivity  is  so  high,  the 
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Fig.  5.  Solution  profiles  of  the  electric  potentials,  electronic  and  ionic  current  fluxes 
in  the  lateral  and  vertical  directions  across  the  MEA  unit  cell  at  lateral  positions  of 
165  |jim,  900  |xm,  and  1500  |xm. 


electric-potential  gradient  within  the  anode  phase  is  very  small. 
Because  the  oxygen-ion  flux  is  dominantly  in  the  direction  nor¬ 
mal  to  the  dense  electrolyte,  the  lateral  electric-potential  gradient 
within  the  electrolyte  phase  (YSZ)  is  also  small.  Upon  transferring 
charge  from  the  oxygen  ion  in  the  YSZ  phase  to  an  electron  in  the 
Ni  phase,  the  Ni  provides  a  low-resistance  path  for  electron  flux 
toward  the  anode  interconnect. 

The  electrode  (Ni)  electric  potential  is  always  lower  than  the 
electrolyte  (YSZ)  potential,  and  the  overpotential  is  always  pos¬ 
itive.  That  is,  the  double  layer  within  the  anode  is  such  that 
A0  =  0Ni  -  @ysz  <  0,  with  the  electric-potential  difference  becom¬ 
ing  less  negative  near  the  dense  electrolyte  interface.  Consequently, 
the  anodic  charge-transfer  rate  is  the  highest  near  the  dense- 
electrolyte  interface.  As  the  magnitude  of  A@  decreases  (i.e.,  the 
YSZ  particles  become  less  positive)  the  rate  of  the  charge-transfer 
reaction  (globally,  H2  (g)  +  O2-  (YSZ)  ^  H20  (g)  +  2e_  (Ni))  increases 
(Fig.  5b).  This  behavior  can  be  understood  qualitatively  in  the  sense 
that  as  the  strength  of  the  double  layer  decreases,  the  barrier 
to  the  transfer  of  negative  charge  into  the  negative  electrode  is 
reduced. 


Cat  lode  current  col 
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_ 

Anode  functional  layer 
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Fig.  6.  Lateral  current  (x  direction,  per  unit  depth)  along  the  width  of  the  MEA. 


Within  the  dense  electrolyte,  there  is  no  need  to  consider 
the  electronic-conducting  electric  potential  because  the  YSZ  is  a 
pure  ionic  conductor.  As  illustrated  in  Fig.  5a,  the  ion-conducting 
electric-potential  profile  within  the  dense  electrolyte  is  linear. 
However,  the  oxygen-ion  flux,  which  is  proportional  to  the  electric- 
potential  gradient  and  the  ion  conductivity,  varies  vertically  within 
the  dense-electrolyte  region  (Fig.  5b). 

Fig.  5c  illustrates  the  electronic  current  density  in  the  lateral 
direction  across  the  width  of  the  MEA  unit  cell.  Because  the  dense 
YSZ  membrane  layer  is  a  pure  ionic  conductor,  the  electronic  cur¬ 
rent  density  vanishes  within  the  dense  YSZ  layer.  Fig.  5c  shows 
that  the  variation  in  lateral  electronic  current  density  through  the 
anode  thickness  is  very  small.  On  the  cathode  side,  however,  the  lat¬ 
eral  electronic  current  density  within  the  LSM  current-collection 
layer  is  much  higher  than  it  is  in  the  LSM-YSZ  layer.  Fig.  5c  also 
shows  that  lateral  current  density  decreases  on  the  cathode  side 
from  the  lateral  position  of  1 65  [xm  toward  1 500  |xm,  and  increases 
correspondingly  at  the  anode  side. 


4. 1.2.  Net  lateral  current 

As  illustrated  in  Fig.  4,  the  distribution  of  the  electrolyte-phase 
electric  potential  and  the  pathlines  of  the  electron  current  fluxes 
within  the  electrodes  indicate  that  the  electron  fluxes  within  the 
anode  and  the  cathode  are  oriented  largely  in  lateral  direction,  and 
the  ion  flux  through  the  dense  electrolyte  is  oriented  nearly  perpen¬ 
dicular  to  the  electrolyte  interfaces  with  the  composite  electrode 
structures.  However,  Fig.  4  does  not  show  quantitatively  the  elec¬ 
tron  and  ion  fluxes,  both  of  which  are  vectors.  The  net  electrical 
currents  flowing  laterally  within  the  anode  and  cathode  structures: 
/a>x  and  Ic<x  are  defined  to  illustrate  the  quantitative  variation  of  the 
electronic  fluxes  as  follows 


fa,x  — 
/c,x  = 


(21) 

(22) 


where  ia<x  and  ic>x  are  the  lateral  components  of  the  electronic  cur¬ 
rent  density  im  (m  =  a  or  c)  and  La  and  Lc  are  thickness  of  the  anode 
and  cathode  layers,  respectively.  The  local  electron  and  ion  current 
densities  can  be  calculated  from  the  electric-potential  gradients  as 


—  —  CTm  V  0  n 


(23) 


Fig.  6  shows  the  net  lateral  electric-current  fluxes  (per  unit 
depth  of  the  cell,  A  cm-1)  within  both  electrodes.  Electric  current 
enters  through  the  cathode  interconnect  and  leaves  through  the 
anode  interconnect.  As  shown  in  Fig.  4,  current  enters  the  cathode 
interconnect  vertically,  then  turns  to  the  lateral  direction.  Thus, 
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Active  MEA  width  (mm) 


Fig.  7.  Profiles  of  net  current  (top  panel),  current  densities  normalized  by  the  overall 
width  of  the  unit  cell  (middle  panel)  and  current  density  normalized  by  the  active 
length  of  the  SIS  unit  cell  (bottom  panel)  as  functions  of  unit-cell  active  width.  The 
results  show  two  cell  voltages,  with  and  without  the  support  structure  in  place.  The 
interconnect  material  has  a  low  electric  conductivity  of  a  =  1  S  cm-1 . 


Active  MEA  width  (mm) 


Fig.  8.  Profiles  of  net  current  (top  panel),  current  densities  normalized  by  the  overall 
width  of  the  unit  cell  (middle  panel)  and  current  density  normalized  by  the  active 
length  of  the  SIS  unit  cell  (bottom  panel)  as  functions  of  unit-cell  active  width.  The 
results  show  two  cell  voltages,  with  and  without  the  support  structure  in  place.  The 
interconnect  material  has  a  high  electric  conductivity  of  a  =  5  x  1 04  S  cm-1 . 


the  net  lateral  electrical  current  increases  near  the  cathode  inter¬ 
connect.  In  the  active  MEA  region,  the  electron  current  within  the 
cathode  is  converted  into  the  ionic  current,  which  flows  through 
the  dense  electrolyte,  and  then  is  converted  into  electron  cur¬ 
rent  within  the  anode.  Because  the  ion  current  through  the  dense 
electrolyte  is  nearly  uniform  (Fig.  5b),  the  lateral  electric-current 
profiles  are  nearly  linear  within  the  active  MEA  regions  of  the 
electrodes  (Fig.  6).  For  this  particular  example,  the  net  current 
developed  by  the  unit  cell  is  about  0.1 6  A  cm-1 .  With  the  cell  oper¬ 
ating  at  Fcen  =  0.75  V,  the  net  power  from  the  cell  is  approximately 
0.1 2  W  cm-1  (per  unit  depth  of  the  cell). 

The  electric  conductivity  of  the  interconnect  structure  can 
greatly  affect  the  cell  performance.  Because  the  interconnect  mate¬ 
rial  must  be  stable  in  both  oxidizing  and  reducing  atmospheres,  the 
materials  choices  are  limited.  One  viable  choice  is  an  electrically 
conducting  doped  perovskite  (e.g.,  strontium-doped  lanthanum 
chromite,  LSC).  However,  because  such  materials  have  relatively 
low  electrical  conductivity  {a  <  1  S  cm-1 ),  there  can  be  a  significant 
voltage  drop  across  the  interconnect.  All  the  electrical  current  pro¬ 
duced  by  the  roughly  1500  pan  width  of  the  active  cell  must  be 
channeled  through  the  relatively  small  interconnect  structure.  This 
current  concentration  can  potentially  cause  local  heating,  or  other 
possibly  deleterious  effects  on  the  cell  materials.  The  design  of  the 
interconnect  regions  and  material  choices  are  important  to  cell  reli¬ 
ability  and  lifetime.  Figs.  7  and  8  show  results  of  a  study  considering 
the  effects  of  increasing  interconnect  electrical  conductivity. 

4.2.  Cell  optimization 


connect  play  important  roles  in  cell  performance.  Consider  all 
geometric  and  physical  properties  to  be  fixed  except  the  unit-cell 
active  width,  which  is  varied  from  0.2  mm  to  1 0  mm.  For  these  com¬ 
parisons  the  fuel  is  taken  to  be  97%  H2  and  3%  H20.  Fig.  7  illustrates 
three  measures  of  cell  performance.  First  is  the  total  current  (per 
unit  depth),  second  is  the  current  density  normalized  by  total  width 
of  the  unit  cell,  and  third  is  the  current  density  normalized  by  the 
active  width  of  the  unit  cell. 

The  total  current  generated  by  the  unit  cell  can  be  evaluated  as 


(24) 


where  “Cl”  is  the  cathode-interconnect  interface  area,  and  ia,y 
is  the  electronic  current  density  in  the  vertical  direction  at  the 
cathode-interconnect  interface,  and  can  be  estimated  based  on 
the  electric  potential  gradient  within  the  interconnect.  The  current 
density  per  active  width  of  the  SIS  unit  cell  can  be  defined  as, 


^active 


I 

^active 


(25) 


where  Wactive  is  the  active  width  of  the  SIS  unit  cell  (Fig.  2).  Simi¬ 
larly,  the  current  density  per  total  width  of  the  unit  cell  is  defined 
as 


itotal_vw 


(26) 


Because  current  in  the  SIS  cell  is  flowing  laterally,  the  lateral 
width  of  the  SIS  unit  cell  and  the  electric  conductivity  of  the  inter¬ 


where  Wtotai  is  the  full  unit-cell  width  (including  the  interconnect 
structure). 
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Fig.  9.  Profiles  of  net  current  (top  panel),  current  densities  normalized  by  the  over¬ 
all  width  of  the  unit  cell  (middle  panel)  and  current  density  normalized  by  the 
active  length  of  the  SIS  unit  cell  (bottom  panel)  as  functions  of  unit-cell  active 
width.  The  results  show  the  effect  of  increasing  electrical  conductivity  of  the  LSM 
current-collection  layer  by  a  factor  of  10.  Two  operating  voltages  are  shown. 

4.2.  I .  Interconnect  conductivity 

Fig.  7  compares  /,  /active,  and  /total  as  functions  of  the  unit¬ 
cell  active  width  at  two  operating  cell  voltages  (Ecen  =  0.5  V  and 
£ceii  =  0.75  V),  with  and  without  the  support  structure  in  place.  In 
this  case  an  interconnect  with  low  electrical  conductivity,  such  as 
LSC,  is  used  (cr  =  l  Son-1).  Fig.  7a  shows  that  the  total  current  I 
(per  unit  depth)  achieves  a  shallow  maximum  at  a  unit-cell  active 
width  of  approximately  7  mm.  When  the  cell  width  is  small  (below 
7  mm  in  this  example),  increasing  the  width  produces  more  total 
current.  Flowever,  as  the  current  increases  the  ohmic  and  concen¬ 
tration  polarization  also  increases.  At  sufficiently  large  unit-cell 
active  width,  the  increased  polarization  causes  the  net  current  to 
begin  decreasing.  Of  course,  the  width  at  which  the  maximum  cur¬ 
rent  is  reached  depends  upon  details  of  the  cell  architecture  and 
materials  properties.  Transport  polarization  through  the  nominal 
support  structure  considered  here  has  a  relatively  small  effect  on 
performance. 

Fig.  7b  illustrates  that  the  current  density  normalized  by  total 
width  of  the  unit  cell  itota This  current  density  is  maximum  at  a 
unit-cell  active  width  of  approximately  1  mm.  Because  the  inter¬ 
connect  structures  are  assumed  to  have  fixed  dimensions,  the 
fraction  of  active  width  compared  to  the  unit-cell  width  decreases 
as  the  unit-cell  width  decreases  (i.e.,  the  interconnect  region  occu¬ 
pies  a  larger  fraction  of  the  unit-cell  width).  At  very  small  unit-cell 
active  width,  itotal  is  small  because  a  large  fraction  of  the  total  unit¬ 
cell  active  width  is  occupied  by  interconnect.  At  large  unit-cell 
active  width,  itotal  decreases  because  ohmic,  activation,  and  concen¬ 
tration  polarizations  increase  as  total  current  increases.  Although 
the  current  increases  with  decreasing  voltage,  the  positions  of  the 
current  maxima  are  affected  only  weakly  by  the  cell  operating 


Fig.  10.  Profiles  of  net  current  (top  panel),  current  densities  normalized  by  the  over¬ 
all  width  of  the  unit  cell  (middle  panel)  and  current  density  normalized  by  the  active 
length  of  the  SIS  unit  cell  (bottom  panel)  as  functions  of  unit-cell  active  width.  A  fuel 
composition  of  15%  H2  and  85%  H20  is  compared  with  a  fuel  composition  of  97%  H2 
and  3%  H20.  Two  operating  voltages  are  shown. 

voltage.  When  high  power  density  for  a  module  is  the  important 
consideration,  itotal  is  a  valuable  measure  of  performance. 

Fig.  7c  illustrates  the  current  density  normalized  by  the  active 
width  /active*  Because  the  interconnect  region  is  excluded,  this  mea¬ 
sure  of  current  density  does  not  have  a  maximum.  As  the  unit-cell 
active  width  increases,  polarizations  associated  with  increasing 
total  current  reduce  the  current  density.  At  large  unit-cell  active 
widths,  the  two  current-density  measures  (itotai  and  iac tive)  con¬ 
verge  toward  the  same  values.  This  is  because  the  interconnect 
region  occupies  a  smaller  relative  fraction  of  the  unit-cell  width. 

Fig.  8  illustrates  the  results  of  a  study  that  is  analogous  to 
that  shown  in  Fig.  7.  In  this  case,  the  interconnect  conductivity  is 
increased  from  cr  =  IS  cm-1  (typical  of  an  electron-conducting  per- 
ovskite)  to  cr  =  5  x  104  S  cm-1  (typical  of  a  noble  metal).  Increasing 
the  interconnect  conductivity  significantly  increases  the  cell  per¬ 
formance  in  the  sense  that  net  current  is  significantly  increased. 
However,  the  position  of  the  maximum  in  total  current  /  is  only 
weakly  affected.  The  current  density  normalized  by  total  width  itotal 
is  increased  to  a  unit-cell  active  width  of  approximately  1.5  mm. 

Although  the  present  simulations  are  based  upon  assuming 
isothermal  operation,  the  local  heating  rates  are  evaluated.  With  an 
interconnect  conductivity  of  o  =  1  S  cm-1  the  local  maximum  heat¬ 
ing  rate  within  the  interconnect  for  the  solution  shown  in  Fig.  4  can 
reach  approximately  2  kW  cm-3.  If  the  interconnect  conductivity  is 
increased  to  cr  =  5  x  104  Scnrr1  the  local  maximum  heating  rate  is 
greatly  reduced  to  approximately  100  Won-3. 

4.2.2.  Cathode  current  collection 

The  electrical  conductivity  of  the  cathode  current-collection 
layer  has  a  large  influence  on  cell  performance.  Although  LSM  is 
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a  relatively  poor  electronic  conductor,  the  LSM  current-collection 
layer  greatly  facilitates  lateral  electron  flux  relative  to  the  LSM-YSZ 
structure  (Fig.  4).  Fig.  9  shows  the  results  of  a  study  that  is  analogous 
to  those  in  Figs.  8  and  7,  where  the  LSM  conductivity  is  increased 
by  a  factor  of  ten.  Although  it  may  not  be  clear  how  to  accomplish 
such  an  increased  conductivity  in  practice,  it  is  clear  the  cell  per¬ 
formance  would  benefit  greatly  by  increasing  the  conductivity  of 
the  cathode  current  collection  layer. 

4.2.3.  Fuel  dilution 

Within  a  full  module  (e.g.,  Fig.  1)  each  unit  cell  is  exposed 
to  a  different  fuel  composition.  Along  the  flow  direction,  the 
fuel  stream  is  degraded  by  consumption  of  the  active  fuels  (e.g., 
H2  and  CO)  and  introduction  of  product  species  (e.g.,  H20  and 
C02).  Fig.  10  illustrates  how  the  optimal  cell  width  varies  as  a 
function  of  fuel  composition.  Although  the  cell  performance  is 
significantly  degraded  as  the  fuel  is  consumed  and  diluted,  the 
optimal  cell  geometry  width  is  only  weakly  affected  by  fuel  com¬ 
position.  This  means  that  the  same  optimal  unit-cell  structure  can 
be  used  for  the  entire  module,  thus  optimizing  the  module  perfor¬ 
mance. 

5.  Summary  and  conclusions 

A  two-dimensional  computational  model  has  been  developed 
to  represent  the  transport  and  chemistry  throughout  a  unit 
cell  of  an  SIS-SOFC  module.  Distributed  electrochemical  charge 
transfer  proceeds  throughout  the  composite  electrode  structures, 
with  Butler-Volmer  charge-transfer  rates  depending  upon  the 
local  gas-phase  composition  and  electric-potential  differences 
between  electrode  and  electrolyte  phases.  Reforming  chemistry 
within  the  anode  structure  is  represented  in  terms  of  a  multi- 
step,  elementary,  reaction  mechanism.  The  model  can  be  applied 
to  evaluate  alternative  cell  architectures  and  materials  proper¬ 
ties. 

Example  problems  are  used  to  illustrate  important  aspects  of 
the  SIS  performance.  The  chemical  and  electrochemical  behav¬ 
ior  is  found  to  be  nearly  one-dimensional  in  the  active  regions 
where  the  dense  electrolyte  is  sandwiched  between  a  cathode  and 
anode  structure.  That  is  to  say,  profiles  of  species  mole  fractions 
and  electric  potentials  depend  primarily  upon  position  normal  to 
the  dense  electrolyte,  but  weakly  upon  the  position  between  the 
cathode  and  anode  interconnectors.  In  the  vicinity  of  the  intercon¬ 
nect  structures,  however,  there  is  strong  two-dimensional  behavior 
as  the  electronic  current  is  channeled  through  the  interconnec¬ 
tion  between  lateral  cells.  Indeed,  these  regions  of  high  electrical 
current  density  deserve  critical  design  consideration,  seeking  to 
mitigate  local  damage  mechanisms. 

The  results  show  that  there  are  optimal  cell  dimensions  that 
maximize  alternative  measures  of  cell  performance.  Increasing  the 
electrical  conductivities  of  the  interconnect  structure  and  the  cath¬ 
ode  current-collection  layer  have  a  significant  beneficial  influence 
upon  cell  performance. 
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